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Abstract. Recently a pair of experiments [1,2] demonstrated a simulation of Abelian anyons in a spin network of 
single photons. The experiments were based on an Abelian discrete gauge theory spin lattice model of Kitaev [3], Here 
we describe how to use hnear optics and single photons to simulate non- Abelian anyons. The scheme makes use of joint 
qutrit-qubit encoding of the spins and the resources required are three pairs of parametric down converted photons and 
14 beam splitters. 
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1. Introduction 

We know that even when the local physics of a system is understood the global behavior can yield surprising 
features that are difficult to predict. Well known examples are critical phenomena and symmetry breaking. 
Topologically ordered media are another example, where one can have a many body system with a strongly 
correlated ground "vacuum" state with more symmetry that the microscopic equations of motion themselves. 
This is manifested in two dimensional electron gases in the quantum Hall effect, where at fractional filling of 
the particle orbitals, the low lying excitations above the ground state are particle-like and predicted to behave like 
anyons [4, 5, 6]. As point-like particles anyons exist only in two dimensions, and unlike Bosons or Fermions which 
have ail phase under exchange, they accumulate an arbitrary phase in the Abelian case, or a matrix valued action 
on the Hilbert space of particle fusion outcomes in the non- Abelian case [7, 8]. 

To date there has been some experimental evidence in support of the existence of Abelian anyons in fractional 
quantum Hall fluids [ /, i u] though there are caveats to interpretation of the data [i 1], and a vigorous experimental 
effort is in place to demonstrate interferometry on non-Abelian anyons. Recently two experiments have reported 
observation of e/4 electron charge in the v — 5/2 filled state [12, 13] which is predicted to have non-Abelian 
quasiparticles. A theoretical interpretation of the data from Ref. [ 1 3] and suggestions for a witness to non-Abelian 
statistics in that experimental setup is given in [14]. There are several proposals for performing non-Abelian 
interferometry in the v = 5/2 filled quantum Hall state [15, 16, 17], and for the general theory of interferometry of 
non-Abelian anyons see Ref s. [18, 19,20]. 

A promising alternative platform to observe anyons is in spin lattices. Here the strategy is to take a two 
dimensional array of spins (i.e. a group of spins whose coupling graph defines a surface) and prepare a highly 
entangled ground vacuum state with respect to some physical theory. Early work by Bais [2 1 ] showed that discrete 
gauge theories in two dimensions could support non-Abelian particles and the full mechanics of braiding in those 
models was given in [22]. Kitaev [3] suggested a spin lattice Hamiltonian that realizes discrete gauge theory on a 
surface with ground states that are invariant under local gauge transformations generated by elements of a discrete 
group. In the case of the simplest Abelian group G = Z2, the microscopic spins are qubits and the model is known 
as the surface code Hamiltonian; 

// = -^A(v)-^B(/) 

/ 

where spins are placed on the edges of a lattice. Here interactions amoung edges meeting at a vertex v [a 
neighborhood denoted star(v)] are given by the vertex operators A (v) = neGstai(v) while the interactions on the 
edges on the boundary of a face / (a neighborhood denoted df) are given by the face operators B{f) = Yleedf'^l- 
The face and vertex terms can be thought of as projections onto gauge invariant states. The ground state has +1 
outcome for measurement of either a face or vertex operator, and a —1 outcome of measuring A(v) or B{f) is 
equivalent to a particle residing on a vertex or face of the lattice; i.e. the particle states are eigenstates of the 
Hamiltonian. 

Several proposals have been made for engineering the surface code Hamiltonian in physical systems, and 
to demonstrate interferometry of Abelian anyons (for a review see [23] and references therein). These proposals 
suggest realizing the Hamiltonian in the perturbative regime corresponding to the gapped phase of native two body 
Hamiltonian on a honeycomb lattice that is realistic [24], and can be designed in the laboratory. With an additional 
magnetic field perturbation, the honeycomb model also supports a phase with non-Abelian anyons and there are 
proposals for observing the statistics in these models [25, 26]. Such simulations require sufficiently large lattices 
to overcome finite size effects [27]. 

Another approach is to not try to build the Hamiltonian at all but rather just simulate the kinematics of the 
anyonic states. That would involve creating the highly entangled ground state, and performing operations on 
the spins for creating, braiding and annihilating excitations. The operations are the same regardless of whether 
the background Hamiltonian is present. Such a digital simulation of the model is not topologically protected since 
there is no gap to creating excitations, but all the same measurement outcomes are obtained provided the operations 
are done quickly relative to decoherence time scales. A promising platform for this kind of digital simulation is to 
use entangled states of single photons. This brings the advantage of easy single qudit manipulation with passive 
optical elements, low decoherence, and well characterised loss channels. Entangling gates between qudits are 
more difficult, but there are well studied pathways for nondeterministically creating such operations. The key 
limitations for optics are a poor single photon sources, limiting experiments to about six photons, and the difficulty 
of building circuits with large depth. Large depth amounts to complex nested interferometers which require tuning 
and stabilisation to a fraction of the wavelength. Recent advances in integrated quantum photonics [28, 29] hold 
the key to increasing the depth of circuits that can be created, and we envisage using just such a platform for 
simulating non-Abelian anyons. 
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It is challenging to build up many body entanglement using single photons and linear optics, yet even the 
small size scale properties are rich. In fact, it is the very virtue of topologically ordered systems that the physics is 
scale invariant. As shown in Fig. 1 a single plaquette of a spin lattice is enough to demonstrate braiding statistics of 
non-Abelian anyons. Of course for such a small system it would be easy to perform a classical simulation to predict 
such behavior. Yet the point would be to exhibit the emergent behavior on small prototypes using tools of coherent 
control that could be extended to scalable systems. Two recent experiments have demonstrated a simulation of 
Abelian statistics in a single plaquette of a spin lattice [1,2] using six and four photons respectively. A simulation 
using NMR control of 4 carbon atoms in a molecule of crotonic acid was also realized [30]. A simulation of non- 
Abelian statistics would mark significant progress for two reasons: first, the effect of non-Abelian statistics is more 
dramatic and easier to discriminate from geometric or dynamical phases, and second, several novel phenomena 
occur such as the existence of topological entropy which could, in principle, be measured in small scale systems. 
The trade off is an increased complexity in the physical operations. Yet for our scheme only three type-I parametric 
down conversion crystals are needed to prepare the initial state, which is the same number used in the Lu et al. [i] 
experiment to measurement AbeUan statistics. 




Figure 1. A simulation of non-Abelian anyonic interferometry using entangled photons (based on a generic protocol in 
[3 1 , 32]). (a) Begin with a surface cellulation where multilevel spins reside on all the links. We could imagine that there 
is a background topological Hamiltonian, Hjo, and the spins are prepared in a highly entangled state corresponding to the 
ground state of //to ■ (b) A single plaquette of the lattice where an operation has been done on one spin, indicated in red, 
which creates localized anionic excitations in //to- The electric charge/anti-charge excitations are drawn as diamonds, 
the flux excitations as squares, (c) Braiding of the flux around one charge, (d) Fusion of the electric charges to measure 
the effect of non-Abelian braiding. An implementation with non-interacting photons would have //to = 0, yet provided 
the entangled vacuum state is prepared and the same operations are done as if //to were present, one obtains the same 
measurement outcomes. 

The manuscript is organized in three stages: the generic lattice model; the mapping to a single plaquette 
system; and the mapping to a physical system. Readers just wishing to see the details of the steps involved in the 
experimental proposal can skip to Sec. 4. In Sec. 2 we introduce the simplest spin lattice model of the discrete 
gauge theories of Kitaev that has non-Abelian anyonic excitations. There are several different types of non-Abelian 
particles arising from this model and in Sec. 3 we describe one of several types of braiding experiments that could 
be done to reveal the non-Abelian statistics. The operations are done using six level particles which can be encoded 
in a qutrit-qubit pair A physical realization of these ideas using linear optical elements and single photons is given 
in Sec. 4. Finally we conclude in Sec. 5 with a discussion of the viability of this technique and possible further 
avenues for experimental inquiry. 

2. Spin lattice model for non-Abelian anyons 

2.1. Non-Abelian discrete gauge theory 

Kitaev's generaUzation of the surface code Hamiltonian is also a sum of vertex and face operators and has localized 
particle-like excitations, but now the ground states are invariant under gauge transformations generated by some 
finite group G = {gj} of choice. We consider the smallest non-Abelian group G = 53, which is the permutation 
group on three objects. To specify the spin lattice we consider a celluation of a two dimensional surface with 
the vertex set V = {v,}, edge set £ = {ej}, and face set J = {//}■ Qudits with d ^\G\=6 levels are placed 

on the edges and physical states are elements of a Hilbert space :H ^ J-C (6)®''^ ' where :}{ (6) = C|0) H h C|5). 

Particles on edges that meet at a vertex v all interact via a vertex operator A (v). Similarly, all particles on edges that 
are on the boundary of a face / interact via B{f ). We pick an orientation for each edge with e = \vj,Vk] denoting 
an edge with arrow pointing from vertex Vj to Vk- The choice of edge orientations is not important as long as a 
consistent convention is used. We assume an orientable complex F and each face / has an orientation consistent 
with it. The Hamiltonian is a sum of operators chosen such that the ground states of Hjq are invariant under local 
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gauge transformations 

7;,(v)= n n (1) 

eyG[v,*] eje[*.y] 

Here Lg{ej),Rg{ej) E U (6), the 6 dimensional unitary group, are the permutation representations of the left and 
right action of multiplication by the group element g E S^on the system particle located at edge ej. For the particle 
states we make the identification \ j) = \gj), where by convention |0) = \go) = \e), with e the identity element. The 
action of left and right group multiplication on the basis states is then L/, |y) = \hgj), and Rh \j) = \gjh) . 
The spin lattice model is: 

//to = -E^M-E^(/) (2) 



where 



{l'k\ne),edf'^k=<'}k\ei,edf 




Figure 2. The triangle of nodes on which to demonstrate the non-Abehan topological action. The white(black) diamonds 
indicate the vertices where the charges(anti-charges) reside, and the aiTows indicate the directions of the edges. The 
dotted ellipses indicate the 6-level qudits on the edges. Each qudit consists of a qutrit, labelled b and indicated by the 
solid circle, and a qubit, labelled a and indicated by the empty circle. The dashed lines indicate the entangled state that 
needs to be created between the three qutiits, and the dotted line indicates the two qubits that must also be entangled. 

In the definition of B{f), the sum is taken over all products of group elements h/^ acting on a counterclockwise 
cycle of edges on the boundary of / such that the accumulated left action is the identity element e G 53 (i.e. 
hiih(_i . . .hjhi = e for the counterclockwise cycle starting at edge e\ and ending at edge e^). The function 
Of{ej) = ±1 according to whether the orientation of the edge is the same as (or opposite to) the face orientation. 
By construction [A(v),A(v')] = [B{f),B{f')] = [A{v),B{f)] = 0. Furthermore, it is straightforward to verify that 
since A(y) is a symmetrized gauge transformation it is a projection, as is B{f). The ground states of Hjq are 
then manifestly gauge invariant states. For a planar surface with boundary, the ground state \GS) of Hjq is unique. 
Excited states can be described by localized particles in the sense that the expectation value (A(v)) = corresponds 
to an excited state with an electric charge located at vertex v while = indicates a magnetic flux at face /. 

These particles have anyonic statistics and we refer the reader to Refs. [31, 32] for a detailed description of the 
particle spectrum of this model and the operations used to create, manipulate and fuse the anyons. 

The specific interferometry experiment we consider is the creation of a fluxless electric charge pair, followed 
by the braiding of a magnetic flux around one charge, and the subsequent fusion of the charges. This entire process 
can be simulated on a single lattice plaquette as summarized in Fig. 1. Pure electric charges are labeled by irreps 
R of the group 53 = {e,c+,c_,fo,fi,f2}, where e is the identity, tj are transpositions, and c± are 3-cycles. The 
non-Abelian electric charges are contained in the two dimensional irrep of 53: 

R2ie)^l2, /?2(rO = crV^^<'^ (c± ) = e±' ^ °\ (4) 
Consider the lattice with a single face and three edges indicated in Fig. 2. The edges are oriented ei = [vi , V2] , 62 = 
[vi, V3],e3 = [v2, V3]. For simplicity of notation we label the edges by the numbers 1,2,3. The ground state is quite 
easy to write down: 

|G5) = i E \k)i®\g)2<»\k-'g)i^W\g)2\k,k-'g),,, (5) 

° S,k€S3 ^ S,k 
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The state of a generic electric charge anti-charge pair located at vertices (vi , V3) in the irrep R (with dimension \R\) 
is given by 

|M«;(vi,V3)) = ^ E ^^{MRR\g)}\g)2 E \k,k-'g)i^i (6) 

withMs an \R\ x \R\ matrix normalized so that Li^j,=i \{MR)a,h\^ — \R\. Under a local gauge transformation T/, at 
vertex vi, the basis states undergo the mapping: Ti,{\'i)\g)2\k,k^^g)i_j, = \hg)2\hk,k^^g)i,j„ and the change in the 
charge state is: 

7/,(vi)|M/j;(vi,V3)) = ^^^,,^sMMRR\g)}\hg)2\hk,k-'g)i3 

= lL,'.kesMMRRHh-'g')}W)2\hk,k-^g)u ,7. 
= lUk'esMMRRHg')m}\8')2\k\k'-'g'),,, ^' 
= \R{h)MR-{vi,V3)). 

where we have used the cyclic property of the trace and the fact that for a unitary representation R^{h^^g') = 
R^ {g')R{h). If instead we had applied the gauge transformation on the other charge we would get 
r/,(v3)|M^;(vi,V3)) = |M^7?(/2^'); (vi, V3)). In fact, these relations are true for any electric charge pair |M^;(v,v')) 
(not necessarily nearest neighbors). 

There is a connection between the action of local gauge transformations and braiding. When a pure magnetic 
flux h E S3 represented as \h) and located at face / is braided in a counterclockwise sense around a charge at 
location v, the action is, 

i?L/j/i)|M«;(y,y')) = |7;(/i)M«; (v, v')>, (8) 

where is the square of the monodromy operator (it represents braiding particle at location / around particle at 
location j in a counterclockwise sense). If a flux with value h is braided around both charges, then the action on 
the state is conjugation: 

^lxl.,\h)\MR;ivy)) = \h)\Rih)MRRih-'y,ivy)), (9) 

There is one state that is invariant under conjugation, the fluxless charge: jl^; (v, v')). 

In fact there is a simple interpretation of this state when the vertices are nearest neighbors, i.e. there is an 
edge e = [v, v']. In that case, the pair of charges is created out of the vacuum by acting on the single spin living on 
the edge e with 

W«(e)|G5) = |lK;(y,v')> (10) 

where 

WR{e) = IgUglxUs) (11) 

with Xk the character of the group element in the irrep R. This operator is in general neither unitary nor Hermitian, 
but for 53 it is Hermitian for all irreps. 

Finally, the amplitude for the process of fusion of the electric charge pair into the vacuum is given by 

FiRih) ^ vac) ^ (l^;(v,v')|/?(/0;(v,v')> = (12) 

The fusion amplitudes can be measured by using an ancillary qubit to perform conditional gauge 
transformations and then measuring the ancilla. For example, say we prepare an ancilla in |+j:)ancj then apply the 
controlled operation 

|0)anc(0|®l + |l)anc(l|®7),(v), (13) 
followed by measurement of the ancilla in the basis | ± v) with outcome m = ztl. The outcome distribution satisfies: 

P(m=l)-P(m=-l)=5R[(l|«|;(v,v')|/?(/i);(v,v'))], 

which is the real part of the fusion amplitude for R{h) 1. Similarly, measuring the ancilla in the basis |±v), 
yields the imaginary part of the fusion amplitude for R{h) vacuum. 

Let's now calculate the expected outcomes for non-Abelian anyonic interferometry. Using Eqs. (4) and (12) 
and the fusion amplitudes to the vacuum are: 

F{R2{e) vac) = 1 

F{R2{tj)^vac) = V; (14) 
F{R2{c±)^vac) = -i. 
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For the measurement proposed here it is also instructive to see how these arise expHcitly in the state overlap. 
The projection operator onto R2 charge is 

WR,=2\e){e\-\c+){c+\-\c^){c4. (15) 
Say we have a fluxless charge pair on nearest neighbor vertices bounding the edge e = [v, v']. Then the overlap is 

{U,;{vy)\R2[hUv,v')) = (l«,;(v,v')|r„|l«,;(v,v'))) 

= {GS\WR,_{e)ThWR^{e)\W) 

= {GS\WR,{e)ThWR,ie)T,^\W) (16) 
= {GS\WR,{e)Wiie)\GS) 
= itr{Wfi,(e)w4>)}. 

Here we have defined the operator Wj{^{e) ~ 2\h)e{h\ — \hc+)e{hc+\ — \hc^)e{hc^\. In the third hne we have used 
the fact that the ground state is invariant under gauge transformations T/, and in the last line we have used the fact 
that the reduced state of the any edge qudit in the ground state is maximally mixed. Evaluating the trace we recover 
the fusion amplitudes in Eq. (14) by considering the state overlap expHcidy. For the given initial state \1r; (v, v')) 
it is also possible to infer a subset of the fusion rules for the theory without using an ancilla at all as we show in 
Sec. 4.3. 



3. Single plaquette scheme 

3.1. Initial state preparation 

An algorithm to prepare the ground state of Hjo in a spin lattice is given in [31]. However, since we only want 
to demonstrate the action of braiding a flux around an electric charge pair, it is sufficient to begin with the state 
I Iff; (^ji''))- Specifically, we want to create the state where the charges are on vertices vi and V3 as in Fig. 2, i.e., 

\lR2;{vi,Vi)) = l[2\e)2(\e,e)n + |fo,fo>i3 + ki,fi>i3 + |f2,f2>i3 + |c+,c_)i3 + \c-,c+)n) 

-|c+)2(|e,C+)l3 + |f0,fl>13 + kl,f2>13 + k2,f0>13 + |c+,e)l3 + |C-,C_)i3) (17) 
-|c_)2(|e,C_)l3 + |fo,f2)l3 + kl,fO>13 + k2,fl)l3 + |c+,C+)i3 + |c_,e)i3)]. 

Here the subscripts indicate the qudit that this is the state for. In Fig. 2 the qudits correspond to the qutrit/qubit 
pairs, and the numbers are shown next to the pair The qutrits and qubits in the pairs are further labeled by "a" for 
the qubits and for the qutrits. 

There are a number of options for how the 6-level qudit can be encoded in the qutrit/qubit pair. It is convenient 
to use the encoding which respects the semidirect product structure of the group ^3 (see Appendix B of Ref. [32]). 
We pick for qudit 1 and 2 the encoding: 

\e) = \l)a\0)b, \t2) = |0)„|2)i, |fo) = |0)„|0)i, |c+) = |1)„|2),, |fi) = |0)„|l)i, |c_) = |1)„|1),,(18) 
and for qudit 3 the encoding: 

\e) = \l)a\0)b, \t2) = \0)a\2)b, ko) = |0)„|0)i, |c+) = |fi) = |0)„|l)fo, \c-) = 1 1)„|2),.(19) 

Using this encoding, the initial state becomes 

\lR,;{vuVi)) = 5|l)2«(|0)l„|0)3„ + |l)l.|l)3«)[2|0)2/.(|0)i,|0)3;, + |l)l^.|l)36 + |2)lfc|2)36) 

-\2)2b{\Q)lb l)3/. + |l)l6|2)3;. + |2)l/,|0)3fe)-|l)2f,(|0)i;.|2)36 + |l)lf,|0)3/. (20) 
+ |2)lfc|l)3/,) . 

This factorises into an entangled pair of qubits and an entangled triple of qutrits, as indicated in Fig. 2. 

Given the initial state preparation, it remains to show how to perform controlled Tc^ and Tt- operations. 

3.2. The controlled Tc^{v\) operations 

The effect of the operation (vi ) is to exchange the basis states of qudits 1 and 2 as 

k) k+), ko) ^ k2), \c+) ^ k-), ki) ^ ko), k-) ^ k), k2) ^ ki>- 

With encoding (18), this corresponds to cyclicly permuting the qutrit states as 

|0)^|2), |1)^|0), |2)^|1). 
This means that, in order to measure the fusion amplitudes and perform the controlled ^^^(vi) operations of 
Eqn. (13), we need to do two controlled permutations — one onlb and one on 2b. 

The operation Tc_ {v\ ) is the inverse, and gives the opposite permutation of the qutrit states 

|0)^|1), |1)^|2), |2)^|0). 
In this case the controlled operation is just the same, except it uses the opposite permutations. 
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3.3. The controlled 7r(,(vi) operation 

If, instead, we want to perform the operation Ttg (vi ), then we need to exchange the states for qudits 1 and 2 as 

k)^ko>, ki)^k+), \t2)^\c-). 

This operation can be achieved by employing a modified CNOT between the qubit and qutrit that only acts on two 
states of the qutrit. These types of modified qubit gates have already demonstrated an advantage in constructing 
multi-qubit gates [33, 34]. First apply such a CNOT with the qubit as control and on states 1 and 2 of the qutrit, 
followed by a NOT gate on the qubit, followed by another CNOT between the qubit and qutrit. Explicitly this 
sequence of operations gives 



k) 


^ |i>« 


0)/.- 


-|l>«|0)i- 


-^|0)„|0)i,- 


- \0)a\0)b ^ 


ko) 


ko) 


= |o>« 


0),,- 


-|0)«|0)i,- 


-^|l)«|0)fc- 


-|l)«|0)i^ 


k) 


ki> 


= |o>« 




-\0)a\l)b- 


-li).-!!)*-- 


-|l)«|2>i^ 


k+ 


h) 


= |0)« 


2)b- 


-\0)a\2)b- 


-^|l)«|2)fo- 


-ll)«ll>i^ 


k- 


k+> 


^ ll>« 


2)h- 


-n)a\l)b- 


-|0>«|l>i,- 


-\0)a\l)b^ 


ki) 


\c-) 


^ ll>« 


1)/.- 


-|l>«|2)i,- 


^\0)a\2)b- 


-|0)„|2)fc^ 





Here the first and third arrows correspond to CNOTs, and the second corresponds to the NOT gate on the qubit. 

To achieve a controlled r,^ (vi ), we can take advantage of the fact that the CNOT applied twice is the identity. 
Therefore, we can apply the CNOT, then apply a CNOT with the ancilla as control and the qubit in the qubit/qutrit 
pair as target, then apply another CNOT between the qubit and qutrit. This needs to be done twice — once for 
qubit/qutrit pair 1 and again for pair 2. A further simplification may be obtained by noting that the measurement 
of the ancilla can take place immediately after the CNOTs with the ancilla as control. The final CNOTs within the 
qubit-qutrit pairs can be omitted, as they do not affect the probabilities of the measurement results. 



3.4. The controlled Tt^ (yi) operation 

Next, the operation Tt^ (vi) exchanges the states for qudits 1 and 2 as 

k>^ki>, k2>^k+>, ko>^k->- 

This operation can be achieved in a similar way to r,„(vi), except with CNOTs acting on states and 1 for the 
qutrit. Explicitly the sequence of operations gives 



k) 


= |1>«|0),- 




-|0)„|1)/.- 


-\0)a\l)b^ 


ki> 


ko) 


= |0)«|0),- 


-^|0)„|0)fe- 


-|1)«|0),- 


-n)a\l)b^ 


k- 


ki) 


= |o)«|i)/,- 


-^|o)«|i)i- 


-nw)b- 


-nU0)b^ 


k) 


k2) 


= |0)«|2)/,- 


-^|0)„|2)i,- 


-|1)"|2)/,- 


-m2)b^ 


k+ 


k+) 


- |1>«|2)/,- 


-^|l)"|2)i- 


-|0)„|2)fc- 


-|0)„|2),^ 


k2> 


k-> 


^ \i)a\l)b- 


-^\l)a\0)b- 


-|0)„|0),- 


-|0)„|0),^ 


ko) 



3.5. The controlled 7;,(vi) operation 

Similarly, the operation TJ^ (vi ) exchanges the states for qudits 1 and 2 as 

k>^k2>, ko>^k+>, ki>^k->- 

Again, the operation follows the procedure above, except with the CNOTs acting on states and 2 for the qutrit. 
Explicitly the sequence of operations gives 



k> 


^ ll>« 


0)/,- 


Hl>«|2)fo- 


-|0)„|2),- 


-|0)„|2),^ 


k2> 


ko) 


= |0>« 


0)/,- 


HO>«|0)fo- 


-|l>«|0)i,- 


A^)a\2)b^ 


k+ 


ki> 


= |o>„ 




HO>„|l)fc- 


-\^)a\l)b- 


-|1>«|1>.^ 


k- 


k2> 


= |o>„ 


2)h- 


H0)„|2)fc- 


-|l>«|2)i,- 


-|1>«|0),^ 


k> 


k+> 


^ |i>« 


2)h- 


Hl>«|0)fo- 


-|0)„|0),- 


- |0)„|0), ^ 


ko) 


k-> 


^ |i>« 


1)/.- 


Hl)«|l>i- 


-|o)«|i)/.- 


-|0)«|l)i^ 


ki> 



Other configurations of the charges turn out to be equivalent to the states above. In the case where both 
charges are on the lower nodes, the encoding on edges 1 and 2 is the same as before, as are the operations that need 
to be performed on them. The only difference is in the encoding of the third edge, which simply corresponds to a 
different interpretation of the experiment. The state obtained if both charges are on the right is again equivalent, 
though this time with different operations. For details see Appendix A. 
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4. Linear optics implementation 

Next we consider how this scheme can be achieved using a linear optics implementation. The core component of 
the linear optical demonstration is the linear optical CNOT of the type demonstrated in Ref. [35]. Such a CNOT 
is based on dual-rail qubits — a single photon in a pair of optical modes. For this gate to function correctly it is 
necessary to ensure that each pair of output modes contains a single photon. This can make it problematic to chain 
these gates together, as it at first appears necessary to perform quantum nondemolition measurements of the photon 
number However, it is still possible provided it can be inferred that a single photon was present in the intermediate 
steps. For example, that is the approach used in Ref. [36]. 

To encode the qutrits we can simply extend the dual-rail encoding to a tri-rail encoding with a single photons 
in three optical modes. This has the advantage that arbitrary unitaries can be performed on these qutrits using 
linear optics [37]. 

Note that for a subset of the fusion rules, we do not need to implement the controlled Tj operations which 
substantially simplifies the experiment. This process is detailed in section 4.3. 



4.1. Initial state preparation 

A possible way to create the three qutrit entangled state required for (vi , V3)), 

{2|0)2i(|0)lfo|0)3i + |l>16|l)3i + |2)l/,|2)3;,)-|2)2/,(|0)ii|l)3/,+ |l)l/.|2)3fc + |2)ii|0)3fc) 

-\l)2b{\0)ib\2hb + \l)ib\Qhb + \2)ib\lhb)}/{2>^), ^ ^ 

is to prepare the initial state |0)2f)|v(/3) \b3b where 

l¥3)ifo3fc = (|0)ii,|0)3/, + |l)ii|l)3fc+|2)ii|2)3i)/%/3 (22) 

is a maximally entangled two qutrit state, and to propagate this state through the circuit in Fig. 3. The phase angles 
(|) and 9 in the circuit are given by 



9 
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7 + \/3' 



2V26 



% 
4' 



Each of the beam splitters depicted acts in a symmetric way, transforming two optical modes described by the 
boson creation operators and b'^ , as i\/R a^ + ^/T~~R and b'^ i^/R b^ + \/T~~R a\ where R is the 
reflectivity indicated. The circuit relies on a final postselection to ensure only one photon is present in each qutrit, 
which happens with a probability of 9/55. The remaining time, invalid qutrit states are produced such as having 
two photons in three modes and these can be postselected out. 

Although the circuit in Fig. 3 looks daunting, each CNOT gate of the type in reference [35] takes 5 beam 
splitter transformations to implement, so the circuit is the same order in complexity as three such gates. In fact, 
the transformation implemented by this circuit is equivalent to a single qutrit transformation on qutrit 2 followed 
by a ternary adder gate between qutrits 2 and 3 \x^x(£iy)). The single qubit transformation would take at 

least two beam splitters and the ternary adder gate can be implemented by four qubit CNOT gates acting between 
pairs of qutrit levels. If such a naive application of qubit gates were used to synthesise the transformation, the 
circuit would already require 22 beam splitters. The circuit in Fig. 3 has been optimised in comparison but it may 
be possible to optimise it further. 

There are potentially a number of ways to prepare the entangled state (22). A fairly direct way, that can 
make use of the same final postselection as the circuit, is to use the output of three type-1 spontaneous parametric 
downconversion (SPDC) crystals in a similar way to Ref. [1]. Each crystal produces an infinite-dimensional 
entangled state between two spatial modes. 



|SPDC) = yT^LX"|««). 

n=0 

For quantum information applications, the state is usually postselected on a single photon being obtained in each 
spatial mode thus producing two photons entangled in energy and momentum. For three independent crystals with 
equal amplitudes for pair production, the state is of the form 

|SPDC)®3 = (l-A.2)3/2 £ A,"i+"2+"3|«i«i,«2n2,n3n3>. 

Here the nj are the photon numbers for the modes produced by crystal j. If we postselect on obtaining a single 
pair of photons, then since the pair may have arisen from any of the three crystals, we obtain the state (22) with 
probabiHty So by using three type-I SPDC sources and an additional photon, and counting three 
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photons in the output, one for each qutrit, we can postselect both the entangled input state (22) and the correct 
operation of the circuit in Fig. 3. 

An alternative method of producing entangled qutrits has been proposed recently [38]. That method has an 
advantage in that it is heralded, but it produces qutrits encoded in the photon-number basis, which makes it more 
difficult to perform subsequent operations. 




Figure 3. The scheme used to produce the desired three qutrit state from a pair of entangled qutrits and an additional 
photon (the state |0)2|i, corresponds to a single photon in mode 2bo). The horizontal lines indicate the different optical 
modes, and the vertical lines indicate symmetric beam splitters connecting modes with reflectivities shown. The open 
circles indicate a relative phase shift in that mode. Note that the beam splitter connecting mode I'o, which is initially 
in a vacuum, and mode Ibi is to introduce loss to balance the circuit. The phase angles (|) and 9, and the beam splitter 
convention, are given in the text. 



4.2. Implementing the controlled Tg operations 

Given our initial state of a pair of electric charge anyons, we can measure the effect of braiding a flux around one 
member of the pair by performing controlled gauge transformations Tg at one vertex (here vertex vi) as described 
in Sec. 2. There are six operators Tg, one for each group element, and recall they involve a product of operations on 
spins that reside on edges that meet at vertex vi . The operator Te is trivial and we first consider how to perform the 
operators T^^ controlled on the state of a qubit ancilla. For these two operators, it is necessary to perform controlled 
swaps on two qutrits in the entangled state, with a qubit as the control. Repeated linear-optical CNOTs with the 
same control qubit would not work, because the linear optical CNOT needs postselection on the final numbers of 
photons. One alternative is to replace the single control qubit with two control qubits, as shown in Fig. 4. The 
CNOTs controlled by the lower rails simply give the required permutation on the first qutrit. The CNOTs with the 
upper rail as control give the inverse permutation. This is followed by the desired permutation without controls, 
with the effect that the desired permutation is performed if the qubit photons are in the lower rails. 

Provided there is a single photon in each qubit or qutrit in the input, and we postselect on single photons at 
the output, then the possibility of photon transfer is eliminated, and the controlled operation is performed correctly. 
For example, if a photon were transferred from the control (qubit 5) to the target (qutrit 1) at the first CNOT, then 
a photon must be transferred from the target (qutrit 1) to the control (qubit 4) at the second to ensure that qutrit 1 
has no more than one photon. However, that would imply that qubit 4 will have more than one photon. 

This convenient procedure will not work if there is the possibility that there are different initial numbers of 
photons at the inputs. For example, if there were initially zero photons in qutrit 1, and two in qutrit 2, then a photon 
could be shifted from qubit 5 to qutrit 1 at the first CNOT, then from qutrit 2 to qubit 5 at the second CNOT. At the 
output, each qubit and qutrit would have a single photon, so this shift would not be detected. If the correct input 
state could be produced without postselection, then this procedure would work. However, for the state preparation 
scheme we have proposed, there is the requirement of postselection of single photons for each output qutrit. 

To achieve the controlled permutation in the case where there is no guarantee of the correct photon numbers 
at the inputs, it would be necessary to use four ancillas in an entangled state. Each CNOT would be applied with 
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Figure 4. The sequence of linear optical CNOTs to peiform the controlled Tc^ operation in the case where there is 
initially a single photon for each qubit and qutrit. The pairs of modes 4 and 5 are the two ancillas, and are initially in an 
entangled state. 



a different ancilla as a control. Because no ancilla is reused, any photon transfer could be detected. Although this 
alternative would work with the states produced by the state preparation scheme, it would be very challenging to 
perform experimentally due to the requirement of generating four entangled photons. 
To achieve the controlled TJ. operation, we can use the following approach. 

(i) Generate the three qutrit entangled state. 

(ii) Apply CNOTs between the qubits and qutrits in pairs 1 and 2. 

(iii) Apply NOTs on the qubits in pairs 1 and 2. 

(iv) Apply CNOTs with an ancilla qubit as control and the qubits in pairs 1 and 2 as targets. 

(v) Measure the ancilla qubit in the ± basis. 

(vi) Perform CNOTs between the qubits and qutrits in pairs 1 and 2 again. 

The final step 6 can be omitted, because it has no effect on the probabilities of the measurement results. In addition, 
the CNOT between 2a and 2/? just simplifies to a NOT, and can be performed deterministic ally. Steps 3 and 4 can 
be combined by making the NOTs controlled by the state (of the ancilla), rather than 1 . The scheme can then be 
achieved as shown in Fig. 5. This example is for 7^^,; the cases of 7}, and Tt^ can be addressed by applying to the 
CNOT and NOT to different pairs of modes in the qutrits. 

There are only three CNOTs required in each case. There is one CNOT between 2a and 2b. Then there is a 
CNOT between the ancilla qubit 4 and la, and between 4 and 2a. These three CNOTs require that there is no photon 
transfer In this case, photon transfer can be ruled out because the ancilla qubits have single photons input, and we 
postselect on single photons output. We also postselect on single photons at the outputs of the qutrits. This ensures 
that there can be no photon transfer, and the CNOTs must have worked. This case is simpler than that for Tc^, 
because only one ancilla qubit is required, rather than four in an entangled state, and one less CNOT is required. 

The anyonic fusion rules given above require that the probabilities for measuring the ancilla qubit in the states 
+ and — are equal. It can immediately be seen that this is what will be obtained, because the CNOT between 4 and 
2a results in these qubits being in an entangled state. As 2a is then discarded without further interaction, this is 
equivalent to just decoherence of the ancilla qubit 4. Then the probabilitites of measuring both + an — are equal 
to 1/2. The difference in probabilities in the case of T^^ is because there is an overlap between the states. 

4.3. Measuring fusion data without controlled operations 

An alternative, simpler, approach to extract information about the non-Abelian statistics is to bypass controlled 
opeations altogether and directly measure operators W^t at edge e, in states of the type \Mr\, where R, R' are 
possibly different irreducible representations. These measurements give nontrivial information about the anyon 
model, namely information about a subset of the fusion rules of the theory. In particular, {Mr\Wj{i\Mr) is always 
zero if R* X R does not contain R'; and by varying M, one can probe the matrices (to be defined below) 

implementing the fusion rules as projectors made of 3 j symbols. 
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qutrit 

















Figure S. The sequence of linear optical CNOTs to perform the controlled T,^ operation. The pair of modes 4 is the 
ancilla, and is initially in a superposition state. The qubit 2a is initially in the state 1 1), so the CNOT just simplifies to a 
NOT operation. 



Let us consider first the case where R = R' = R2, the two-dimensional irrep of S3. In this case, we do 
not need to perform controlled operations T,, or Tc+, and just need the corresponding unconditional operations. 
To see how this works we compute the outcome of the measurement of in the gauge transformed state 
7)Xv)|l«,;(v,v')>-|/?2W;(v,v')>, 



{R2{hy,{vy)\WR,{e)\R2{h)-{vy)) 



(l«,;(v,v')|7;>«2(e)7),(v)|l«2;(v,v'))) 
{GS\WR,{e)Wli;'{e)WR,{e)\GS) 

itr{(H'^,(.))X:'W}- 



We find 



{R2ihy,ivy)\WR,ie)\R2ih);ivy)) = 



1 = , 

1 

2 



— 4 /i = c± 



(23) 



(24) 



h = tj'ij 

This reproduces the fusion amplitudes computed before, but this is just an accident. In fact, such an expectation 
value probes the fusion rules of the theory: 

\R\ \R'\ 



{Rih)\WR,\R{h)) = l.Qai;jK>,mcie{h) , 

a,b,d,e= \ c=\ 



(25) 



where are the projectors onto the vacuum fusion channel for three irreducible representations, that is, 



g a 

where the gr[^'''^'^'^*'']-" are an orthonormal basis for the +1 eigenspace of gl^ ■'^"■'^l; they are more familiar as 3 j 
symbols in the case of angular momentum. 

Measurements for {R2{h); {v,v')\WR2{e)\R2{h); {v,v')) can easily be performed deterministically. T^j^ isjusta 
permutation of the three modes on the qutrit. The operation Tfg gives the change in the encoded states 

|l>„|0)/,^|0)„|0)fo, |0)„|l)fo^|l)„|2),, |0)„|2)fo^|l)„|l),. 

This can be achieved by a NOT gate on the qubit, and swapping states 1 and 2 for the qutrit. Then operation r,(, 
is the same, except swapping states and 1 for the qutrit, and Ttg is the same except swapping states and 2. 
These are all operations that can be performed deterministically. Then the operator Wr2 can just be measured by 
measuring the qubits and qutrits in their computational basis. 

In particular, with the operation T^^ the unnormalised state becomes 

|l>2«(|0>l„|0)3„+|l)l„|l)3„)[2|2>2i,(|2>lfc|0)3fo+|0)i,|l)3i+|l)l/,|2)3i,) 
-|l)2fo(|2)l/,|l)3i,+ |0)ifo|2)3i,+ |l)lfo|0)3i,)-|0)2i,(|2)ii|2)3fo+|0)l,,|0)3i + |l)l/,|l)3fo)]. 



(26) 
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Then the probabilities of measuring e, c+ and c_ are 1/6, 2/3 and 1 /6 respectively, which gives (Wrj) = —1/2- 
With the operation Tt^ the unnormalised state becomes 

|0)2«(|l)u,|0)3„+|0)l„|l)3„)[2|0)2fo(|0)li,|0)3/,+ |2)ii,|l)3fe + |l)l/,|2)3i,) 

-|l)2/,(|0)l/,|l)3i,+ |2)li|2)3i,+ |l)lfo|0)3i,)-|0)2i,(|0)ii,|2)3i+|2)ii,|0)3i+|l)l/,|l)3fo)]. ^ ' 

In this case, the state \Q)2a means that the probabilities of e, c+ and c_ are all zero, so (Wr,) = 0. 

On the other hand, for the original state the probabilities of e, c+ and c_ are 2/3, 1/6 and 1 /6, respectively, 
giving (Wrt) = 1- This approach gives us a method to test the system deterministically, once the correct initial 
entangled state is prepared. 

5. Conclusions 

In summary we have proposed a physical implementation of non-Abelian interferometry using entangled multi- 
photon states. The realization is via a minimal construction of a non-Abelian discrete gauge theory in a single 
plaquette of a spin lattice. The entire construction requires three parametric down conversion crystals to generate 
entangled photon pairs which are then processed using beam splitters and phase shifters to simulate the braiding 
and fusion of non-Abelian charges. By choosing an initial state which is already an excited state of the model 
with a pair of vacuum electric charge pairs, and by judicious choice of witness to the fusion data we have found a 
substantially simplified protocol for the interferometry. 

Given recent advances in integrated photonic devices our protocol should be within experimental reach in the 
relatively near future. Possible extensions of this model include the simulation of thermal states of topologically 
ordered media by introducing mixedness using ancillary degrees of freedom (see e.g. [39]). Such measurements 
could then be used to explore phenomena like thermal fragility of topological entanglement [40]. 
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Appendix A. Other charge configurations 

Alternative configurations of the charges turn out to be equivalent. Consider the case that the charges are on the 
lower two nodes. The state in this case is 

\Ir2,{vIjV2)) = ^[2\e)i(\e,e)23 + |fO,fo)23 + kl,ri)23 + |r2,f2)23 + |c+,C+)23 + |C-,C_)23) 

-|c+)i (|e,C_)23 + |f0,fl>23 + kl,f2>23 + k2,fo)23 + |c+,e)23 + |C-,C+)23) (A.l) 
-|c_)i (|e,C+)23 + ko,f2>23 + kl,fo)23 + k2,fl)23 + k+,C_)23 + |c_,e)23)]- 

In this case, we can use the encoding (18) on all three qudits, and the state becomes the same as (20), except with 
the labels 1 and 2 interchanged. 

Another case is where both charges are on the right, and the state becomes 

|1«2;(^2,V3)) = 5[2|e)3(|e,e)l2 + |fo,fo>i2 + kl7fl>12 + 1^2, ^2) 12 + |c+,C+)l2 + \C-,C-)l2) 

-\c+)J\e,c+)n + ko,fi>i2 + ki,f2>i2 + k2,fo)i2 + k-,e)i2 + \c+,c^)n) (A.2) 

-k->3(k,C-)l2 + k0)f2)l2 + kl)fo)l2 + k2,fl>12 + |C-,C+)i2 + |c+,e)l2)]. 

Here we would use the encoding (18) on qudit 3, and the encoding (19) on qubits 1 and 2. Then we would again 
get the state (20), except this time with the labels 2 and 3 interchanged. 

The encodings are different on the nodes we wish to perform the operations on, but the operations are different, 
and the scheme is again isomorphic. To achieve the operations Tg{v2) we need to apply operations to qudits 1 and 
3. The encoding on 3 is now the same as the encoding on 1 and 2 was previously, and we again are performing L 
operations on this qudit, so the analysis is identical to that above. For qudit 1 , we now need to perform R operations, 
because this edge is directed inwards towards node 2. The effect of the operation T^^ (V2) is to exchange the basis 
states of qudit 1 as 

k> ^ k->, ko> ^ k2>, k-> k+>, ki> ko>, k+> ^ k>, k2> ^ ki>- 

In this case, because we are now using the encoding (19) on qudit 1, this corresponds to permuting the basis states 
as 

|0)^|2), |1)^|0), |2)^|1). 
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This is identical to that obtained previously. The requirement to perform the R operation, together with the different 
encoding, mean that the end result is identical. 

The results are similar for the other cases. Tc_{v2) is just the inverse, and the 7r,(vi) are also the same as 
before. For example, to perform the operation TJq (vi ), then we need to exchange the states for qudit 1 as 

k)^ko), ki>^|c-), k2>^|c+). 

Again, the roles of c_ and c+ are exchanged, but the encoding also exchanges the roles of c_ and c+, so the 
operations required on the encoded states are identical. 

Recall also that in this case, the roles of qudits 2 and 3 are exchanged for the state, so operations on qudits 
1 and 3 here are equivalent to operations on qudits 1 and 2 previously. As a result of this, in each of the three 
possibilities for the charges, both the encoded states and the operations on these encoded states are the same. 
Hence any experiment can be interpreted in each of three ways: as with the charges on the nodes shown in Fig. 2, 
as with the charges on the bottom two nodes, or with the charges on the right two nodes. 
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